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Abstract. We present a novel algorithm for deciding whether a given planar curve is an 
image of a given spatial curve, obtained by a central or a parallel projection with unknown 
parameters. The motivation comes from the problem of establishing a correspondence be- 
tween an object and an image, taken by a camera with unknown position and parameters. 
A straightforward approach to this problem consists of setting up a system of conditions 
on the projection parameters and then checking whether or not this system has a solution. 
The computational advantage of the algorithm presented here, in comparison to algorithms 
based on the straightforward approach, lies in a significant reduction of a number of real 
parameters that need to be eliminated in order to establish existence or non-existence of 
a projection that maps a given spatial curve to a given planar curve. Our algorithm is 
based on projection criteria that reduce the projection problem to a certain modification 
of the equivalence problem of planar curves under affine and projective transformations. 
To solve the latter problem we make an algebraic adaptation of signature construction that 
has been used to solve the equivalence problems for smooth curves. We introduce a notion 
of a classifying set of rational differential invariants and produce explicit formulas for such 
invariants for the actions of the projective and the affine groups on the plane. 

Key words: central and parallel projections; finite and affine cameras; camera decom- 
position; curves; classifying differential invariants; projective and affine transformations; 
signatures; machine vision 



1 Introduction 

Identifying an object in three-dimensional space with its planar image is a fundamental problem 
in computer vision. In particular, given a database of images (medical images, aerial pho- 
tographs, human photographs), one would like to have an algorithm to match a given object 
in 3D with an image in the database, even though a position of the camera and its parameters 
may be unknown. Since the defining features of many objects can be represented by curves, 
obtaining a solution for the identification problem for curves is essential. 

A central projection from M 3 to M 2 models a pinhole camera pictured in Fig. 1. It is described 
by a linear fractional transformation 
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Figure 1. A pinhole camera, http://en.wikipedia.0rg/wiki/File:Pinh0le-camera.png. 

where (z±, z<i,z-$) denote coordinates in M 3 , (x, y) denote coordinates in M 2 and p^, i = 1, . . . , 3, 
j = 1, ... ,4, are real parameters of the projection, such that the left 3x3 submatrix of 3 x 4 
matrix P = (pij) has a non-zero determinant. Parameters represent the freedom to choose 
the center of the projection, the position of the image plane and (in general, non-orthogonal) 
coordinate system on the image plane 1 . In the case when the distance between a camera and 
an object is significantly greater than the object depth, a parallel projection provides a good 
camera model. A parallel projection has 8 parameters and can be described by a 3 x 4 matrix 
of rank 3, whose last row is (0, 0, 0, 1). We review various camera models and related geometry 
in Section 2 (see also [14, 20]). In most general terms, the object-image correspondence problem, 
or the projection problem, as we will call it from now on, can be formulated as follows: 

Problem 1. Given a subset Z of IR 3 and a subset X of M 2 ; determine whether there exists 
a projection P : M 3 > R 2 , such that X = P(Z) ? 2 

A straightforward approach to this problem consists of setting up a system of conditions on 
the projection parameters and then checking whether or not this system has a solution. In the 
case when Z and X are finite lists of points, a solution based on the straightforward approach 
can be found in [20]. For curves and surfaces under central projections, this approach is taken 
in [15]. However, internal parameters of the camera are considered to be known in that paper 
and, therefore, there are only 6 camera parameters in that study versus 12 considered here. The 
method presented in [15] also uses an additional assumption that a planar curve X C R 2 has at 
least two points, whose tangent lines coincide. An alternative approach to the problem in the 
case when Z and X are finite lists of points under parallel projections was presented in [1, 2]. 
In these articles, the authors establish polynomial relationships that have to be satisfied by 
coordinates of the points in the sets Z and X in order for a projection to exists. 

Our approach to the projection problem for curves is somewhere in between the direct ap- 
proach and the implicit approach. We exploit the relationship between the projection problem 
and equivalence problem under group-actions to find the conditions that need to be satisfied by 
the object, the image and the center of the projection 3 . In comparison with the straightforward 
approach, our solution leads to a significant reduction of the number of parameters that have to 
be eliminated in order to solve Problem 1 for curves. 

All of the theoretical results of this paper are valid for arbitrary irreducible algebraic curves 
(rational and non-rational), but the algorithms are presented for rational algebraic curves, 

It is clear from (1) that multiplication of P by a non-zero constant does not change the projection map. 
Therefore, we can identify P with a point of the projective space PR 11 , rather than a point in R 12 . However, 
since we do not know which of the parameters are non-zero, in computations we have to keep all 12 parameters. 

2 We borrow the notation $ : V > W for a rational map from V to W from the algebraic geometry literature 

(see for instance, [12]) in order to emphasize that the map is defined almost everywhere on V . We use the same 

letter P to denote the 3 x 4 matrix P = (pij) and the map P: R 3 > R 2 defined by (1). Hence P{Z) = {(x,y) G 

R 2 I 3 z G Z such that (1) holds}. 

3 In the case of parallel projection, when the center is at infinity, the conditions are on the direction of the 
projection. 
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i.e. Z = {r(s) | s in the domain of T} and X = {"f(t) \ t in the domain of 7} for rational maps 

T: M > M 3 and 7: K > K 2 . A bar above a set denotes the Zariski closure of the set 4 . 

Throughout the paper, we assume that Z is not a straight line (and, therefore, its image under 
any projection is a one-dimensional constructible set). Since, in general, P{Z) is not Zariski 
closed we must relax the projection condition to P{Z) = X . Under those conditions, Problem 1, 
for central projections, can be reformulated as the following real quantifier elimination problem: 

Reformulation 1 (straightforward approach). Given two rational maps T and 7, determine 
the truth of the statement: 

]PeWcK 3x4 Vs in the domain of T 3teR P(T(s)) = j(t), 

where U is the open subset of the set of 3 x 4 matrices defined by the condition that the left 3x3 
minor is nonzero 5 . 

Real quantifier elimination problems are algorithmicaily solvable [30]. A survey of subse- 
quent developments in this area can be found, for instance, in [22] and [11]. Due to their high 
computational complexity (at least exponential) on the number of quantified parameters, it is 
crucial to reduce the number of quantified parameters. The main contribution of this paper is 
to provide another formulation of the problem which involves significantly smaller number of 
quantified parameters. 

We first begin by reducing the projection problem to the problem of deciding whether the 
given planar curve X is equivalent to a curve in a certain family of planar curves under an 
action of the projective group in the case of central projections, and under the action of the 
affine group in the case of parallel projections. The family of curves depends on 3 parameters 
in the case of central projections, and on 2 parameters in the case of parallel projections. 

Then we solve these group-equivalence problems by an adaptation of differential signature 
construction developed in [9] for solving local equivalence problems for smooth curves. We give 
an algebraic formulation of the signature construction and show that it leads to a solution of 
global equivalence problems for algebraic curves. For this purpose, we introduce a notion of 
a classifying set of rational differential invariants and obtain such sets of invariants for the 
actions of the projective and affine groups on the plane. Following this method for the case 
of central projections, when Z and X are rational algebraic curves, we define two rational 
signature maps S\x- — > K 2 and S\z- M 4 — — > M 2 . Construction of these signature maps 
requires only differentiation and arithmetic operations and is computationally trivial. Then 
Problem 1 becomes equivalent to 

Reformulation 2 (signature approach). Given two rational maps S\x and S\z, determine the 
truth of the statement: 

3 c £ W C I 3 V s in the domain of Sz (c, s) 3 1 G E Sz (c, s) = Sx it) , 

where U is a certain Zariski open subset 0/M 3 . 

Note that Reformulations 1 and 2 have similar structure, but the former requires elimination 
of 14 parameters (pn, . . . ,P3i,s,t), while the latter requires elimination of only 5 parameters 
(ci, C2, C3, s, t). The case of parallel projection is treated in the similar manner and leads to the 
reduction of the number of real parameters that need to be eliminated from 10 to 4. 

4 Recall that a set W C R n is Zariski closed if it equals to the zero set of a system of polynomials in n variables. 
The complement of a Zariski closed set is called Zariski open. A Zariski open set is dense in R n . A Zariski 
closure W of a set W is the smallest (with respect to inclusions) Zariski closed set containing W. 

5 Note that, in Reformulation 1, we decide whether P(Z) C X, which appears to be weaker than P(Z) = X. 
However, they are actually equivalent. Since we assumed that Z is not a line, the set P(Z) is one-dimensional. 
Since X is rational algebraic curve, it is irreducible. Hence P(Z) C X <^=> P{Z) = X. 
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Although the relation between projections and group actions is known, our literature search 
did not yield algorithms that exploit this relationship to solve the projection problem for curves 
in the generic setting of cameras with unknown internal and external parameters. The goal of 
the paper is to introduce such algorithms. The significant reduction of the number of parameters 
in the quantifier elimination problem is the main advantage of such algorithms. 

A preliminary report on this project appeared in the conference proceedings [8] . The current 
paper is significantly more comprehensive and rigorous, and also includes proofs omitted in [8]. 
Although the development of efficient implementation lies outside of the scope of this paper, we 
made a preliminary implementation of an algorithm based on signature construction presented 
here and an algorithm based on the straightforward approach over complex numbers. The Maple 
code and the experiments are posted on the internet [7]. The existence of a projection over com- 
plex numbers provides necessary but not sufficient condition for existence of a real projection. 

The paper is structured as follows. In Section 2, we review the basic facts about projections 
and cameras. In Section 3, we prove projection criteria that reduce the central and the parallel 
projection problems to a certain modification of the projective and the affine group-equivalence 
problems for planar curves. This criteria are straightforward consequences of known camera de- 
compositions [20]. In Section 4, we define the notion of a classifying set of rational differential 
invariants and present a solution of the global group-equivalence problem for planar algebraic 
curves based on these invariants. This is an algebraic reformulation of a solution of local group- 
equivalence problem for smooth curves [9]. In Section 5, combining the ideas from the previous 
two sections, we present and prove an algorithm for solving the projection problem for rational 
algebraic curves and give examples. In Section 6, we discuss possible adaptations of this algo- 
rithm to solve projection problem for non-rational algebraic curves and for finite lists of points. 
We discuss the subtle difference between the discrete (with finitely many points) and the con- 
tinuous projection problems, showing that the solution for the discrete problem does not provide 
an immediate solution to the projection problem for the curves represented by samples of points. 
This leads us into the discussion of challenges that arise in application of our algorithms to real- 
life images, given by discrete pixels, and of ideas for overcoming these challenges. In Appendix A, 
we give explicit formulae for affine and projective classifying sets of rational invariants. 

2 Projections and cameras 

We embed R n into projective space PR n and use homogeneous coordinates on PR n to express 
the map (1) by matrix multiplication. 

Notation 1. Square brackets around matrices (and, in particular, vectors) will be used to 
denote an equivalence class with respect to multiplication of a matrix by a nonzero scalar. 
Multiplication of equivalence classes of matrices A and B of appropriate sizes is well-defined by 
[A][B]:=[AB]. 

With this notation, a point (x, y) G R 2 corresponds to a point [x, y, 1] = [Ax, Xy, A] G PR 2 for 
all A / 0, and a point {z\, z 2 , zs) G R 3 corresponds to [z\, z 2 , z%, 1] G PR 3 . We will refer to the 
points in PR™ whose last homogeneous coordinate is zero as points at infinity. In homogeneous 
coordinates projection (1) is a map [P] : PR 3 — > PR 2 given by 

[x,y,l] T = [P][z 1 ,z 2 ,z 3 ,l] T , 

where P is 3 x 4 matrix of rank 3 and superscript T denotes transposition. Matrix P has a 1- 
dimensional kernel. Therefore, there exists a point [z®, z\, z®, z®] G PR 3 whose image under the 
projection is undefined (recall that [0,0,0] is not a point in PR 2 ). Geometrically, this point is 
the center of the projection. 
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In computer science literature (e.g. [20]), a camera is called finite if its center is not at 
infinity. A finite camera is modeled by a matrix P, whose left 3x3 submatrix is non-singular. 
Geometrically, finite cameras correspond to central projections from M 3 to a plane. On the 
contrary, an infinite camera has its center at an infinite point of PR 3 . An infinite camera is 
modeled by a matrix P whose left 3x3 submatrix is singular. An infinite camera is called affine 
if the preimage of the line at infinity in PR 2 is the plane at infinity in PR 3 . An affine camera 
is modeled by a matrix P whose last row is (0, 0, 0, 1). In this case map (1) becomes 

X = PllZi + pi 2 Z 2 + P13Z3 + PU, V = P2\Z\ + P22Z2 + P23Z3 + P24- 

Geometrically, affine cameras correspond to parallel projections from M 3 to a plane 6 . Eight 
degrees of freedom reflect a choice of the direction of a projection, a position of the image plane 
and a choice of linear system of coordinates on the image plane. In fact, by allowing the freedom 
to choose a non-orthogonal coordinate system on the image plane, we may always assume that 
we project on one of the coordinate planes. 

Definition 1. A set of equivalence classes [P], where P is a 3 x 4 matrix whose left 3x3 
submatrix is non-singular, is called the set of central projections and is denoted CV. 

A set of equivalence classes [P], where P has rank 3 and its last row is (0,0,0, A), A ^ 0, is 
called the set of parallel projections and is denoted VV. 

Equation (1) determines a central projection when [P] G CV and it determines a parallel 
projection when [P] G VV. Sets CV and VV are disjoint. Projections that are not included in 
these two classes correspond to infinite, non-affine cameras. These are not frequently used in 
computer vision and are not considered in this paper. 

3 Reduction to the group-equivalence problem 

Definition 2. We say that a curve ZcK 3 projects to X C M 2 if there exists a 3 x 4 matrix P 
of rank 3 such that X = P(Z), where 

P(Z) = {{x,y) G M 2 I 3z G Z such that (1) holds}. 

Recall that for every algebraic curve X C W l there exists a unique projective algebraic 
curve [X] C PM n such that [X] is the smallest projective variety containing X (see [17]). 
It is not difficult to check that X = P{Z) is equivalent to [X] = [P][Z], where [-P][i?] = 
{[P][z]\[z]e[Z], [P][z]^[0]}. 

Definition 3. The projective group 7 VQC{n + 1) is a quotient of the general linear group 
QC{n + 1), consisting of (n + 1) x (n + 1) non-singular matrices, by a 1-dimensional abelian 
subgroup XI, where A / £ K and / is the identity matrix. Elements of VQC(n + 1) are 
equivalence classes [B] = [XB], where A 7^ and B G QC{n + 1). 

The affine group A{n) is a subgroup of VQC{n + l) whose elements [B] have a representative 
B G QC(n + 1) with the last row equal to (0, ... , 0, 1). 

The equi-affine group SA{n) is a subgroup of A{n) whose elements [B] have a representative 
B G QC{n + 1) with determinant 1 and the last row equal to (0, . . . , 0, 1). 

In homogeneous coordinates, the standard action of the projective group VQC{n + 1) on PR n 
is defined by multiplication 

[zi, . . . , z n , z ] T -)■ [B] [zi, . . . , z n , z ] T . (2) 

6 Parallel projections are also called generalized weak perspective projections [1, 2]. 

7 We will occasionally include a field in the group-notation, e.g. VGC(n,C) or VGC(n, R). If the field is not 
indicated we assume that the group is defined over R. 
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The action (2) induces linear-fractional action of VQC(n + 1) on M n . 8 The restriction of (2) 
to A(n) induces an action on W n consisting of compositions of linear transformations and trans- 
lations. 

Definition 4. We say that two curves X\ C W 1 and X 2 C W 1 are VQC(n + l)-equivalent if there 

exists [A] £ VQC{n + 1), such that [X 2 ] = {[A][p] | [p] £ [Xi]}. We then write X 2 = A ■ X 1 or 

[X2] = [A][X\]. If [A] £ G, where G is a subgroup of VQC[n + 1), we say that X\ and X 2 are 

G-equivalent and write X\ = X%. 

G 

Before stating the projection criteria, we make the following simple, but important observa- 
tions. 

Proposition 1. 

(i) If Z C R 3 projects to X C R 2 by a parallel projection, then any curve that is .4(3)- 
equivalent to Z projects to any curve that is A(2) -equivalent to X by a parallel projection. 
In other words, parallel projections are defined on affine equivalence classes of curves. 

(ii) //2 C I 3 projects to X C R 2 by a central projection then any curve in R 3 that is .4(3)- 
equivalent to Z projects to any curve on R 2 that is VGC(3) -equivalent to X by a central 
projection. 



Proof, (i) Assume that there exists a parallel projection [P] £ VV such that [X] 
Then for all (A, B) G A{2) x ,4(3) we have 



[P\[Z]. 



[A] [X] = [A] [P] [B-t] [B] [Z] = [A] [P] [B-t] [B] \Z\. 

Since [A] [P] [B^ 1 ] £ VV, curve B • Z projects to A ■ X. (ii) is proved similarly. ■ 

Remark 1. It is known that if X\ and X 2 are images of a curve Z under two central projections 
with the same center, then X\ and X 2 are 'P(/£(3)-equivalent, but if the centers of the projections 
are not the same this is no longer true (see Example 3). Similarly, images of Z under various 
parallel projections may not be „4(2)-equivalent. 

Theorem 1 (central projection criterion). A curve Z C M 3 projects to a curve X C M 2 by 
a central projection if and only if there exist c\,c 2 ,c^ £ M such that X is VQ X(3) -equivalent to 
a planar curve 



Z, 



'Cl,C2,C 3 — 



' Z\ + Cl z 2 + c 2 
K z 3 + c 3 ' z 3 + c 3 



[zi,Z 2 ,Z3) £ Z 



(3) 



Proof. (=^>) Assume there exists a central projection [P] such that X = P(Z). Then P is a 3 x 4 
matrix, whose left 3x3 submatrix is non-singular. Therefore there exist c\, c 2 , C3 £ M such that 
p*4 = cip*i + c 2 p* 2 + C3f>*3, where p*j denotes the j-th column of the matrix P. We observe that 

[A][P$][B] 



[P], 



(4) 



where A is the left 3x3 submatrix of P, 



r c ■- 



'I N 
10 
,0 10, 



and 



B :-- 



(I ci\ 

1 c 2 

1 c 3 

\0 1 ) 



(5) 



8 Linear-fractional action of VQLin + 1) on R n is an example of a rational action of an algebraic group on an 
algebraic variety. General definition of a rational action can be found in [23, Definition 2.1]. 
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Note that [A] belongs to PQL{3). Since 

[Pg][B][z 1} Z2,Z3, 1] T = [Zl + C!, Z 2 + C 2 , Z 3 + C 3 ] T , 

then [X] = [A][Z C1 ,c2,c 3 ]) where Z Cl:C2fi3 is defined by (3). 

(<=) To prove the converse direction we assume that there exists [A] G PQC(3) and 01,02,03 G 
K such that [Af] = [-A][>Z cl)C2jC3 ], where Z cljC2jC3 is defined by (3). A direct computation shows 
that Z is projected to X by the central projection [P] = [A] [P®] [B], where B and [P®] are given 
by (5). ■ 

We note that the map 

Zl + Cl z 2 + c 2 

& = : , V= : 6 

Z3 + c 3 Z3 + c 3 

is a projection centered (— c\, —C2, — 03) to the plane 213 = 1 with coordinates on the image plane 
induced from M 3 , namely, x = z\ and y = Z2- We call (6) the canonical projection centered at 
(—ci, —C2, — C3). It follows from decomposition (4) that any central projection is a composition 
of a translation in M 3 (corresponding to translation of the camera center to the origin), the 
canonical projection P® centered at the origin, and a projective transformation on the image 
plane. 

Remark 2 (CP is a homogeneous space). It is easy to check that the map 

* : {VGC{3) x .4.(3)) xCP^CP 
defined by 

^(([A},[B}),[P]) = [A][P}[B- 1 ] (7) 

for [P] G CP and ([A], [B]) G -P££(3) x 4.(3) is an action of the product group PQC{3) x A{3) 
on the set of central projections CP. Decomposition (4) shows that this action is transitive. The 
stabilizer of the canonical projection P® centered at the origin is a 9-dimensional group 



H u c = i[[A] 1 °' 



where A G QC(3). 



1 

The set of central projections CP is, therefore, diffeomorphic to the homogeneous space 
PGC{3) x A(3)/H$. 

Theorem 2 (parallel projection criterion). A curve Z C M 3 projects to a curve X C M 2 
by a parallel projection if and only if there exist c\,C2 G M and an ordered triplet (i,j,k) G 
{(1, 2, 3), (1, 3, 2), (2, 3, 1)} such that X is 4,(2) -equivalent to 

Zclc2 = + c i^> z 3 + c 2 z fc) I ( z 1^2,z 3 ) G Z}. (8) 



Proof. (=>) Assume there exists a parallel projection [P] such that = P(Z). Then [P] can 
be represented by a matrix 

'v\\ P12 P13 Pu' 
P21 P22 P23 P2A 
1 
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of rank 3. Therefore there exist 1 < i < j < 3 such that the rank of the submatrix 



Pit Plj\ 

P2i P2j) 

is 2. Then for 1 < k < 3, such that k ^ i and k ^ j, there exist ci,c 2 G M, such that 

Pli Plj P14 N 



Plfe 
P2k 



Cl 



+ c 2 



Plj 

P2j 



matrix whose columns are vectors b* 
6*4 = (0, 0, 0, 1) T . We observe that 



We define A := 

(1,0,0,0) T , 



P2i P2j P24 
1 



[P] = [A][P%][B], 



where P v :- 



Since [P%][B][Z] = then [X] 



1 N 
10 
,0 1, 



and define -B to be the 
(0,1,0,0) T , := (ci,c 2 ,l,0) T , 

(9) 



[^][Zcf; C2 ]. Observe that [A] £ A{2) and the direct 



statement is proved. 

(<=) To prove the converse direction we assume that there exist [A] £ .4(2), two real num- 
bers ci and c 2 , and a triplet of indices such that (i, j, k) £ {(1, 2, 3), (1, 3, 2), (2, 3, 1)}, such that 



[X] = [A]{Z l c ^' C2 ], where a planar curve Zc'i,'c 2 is given by (8). Let B be a matrix defined in the 
first part of the proof. A direct computation shows that Z is projected to X by the parallel 
projection [P] = [A][F$\[B]. ■ 

Remark 3 (VV is a homogeneous space). The map ^ : (A(2) x -4(3)) x TV — > VV defined 
by (7) for [P] £ VV and ([A], [B]) £ A(2) x -4(3) is an action of the product group -4(2) x -4(3) 
on the set of parallel projections VV. Decomposition (9) shows that this action is transitive. 
The stabilizer of the orthogonal projection Pp is a 10-dimensional group 



f / 



mn mi2 a\ 
77121 m 22 02 
1 



"in Tn\2 a\ 

77Z21 ™-22 a 2 

m 31 m 32 m 33 a 3 

1 



\ 



/ 



m 33 (miim 2 2 - mi 2 m 2 i) / 



The set of central projections "P"P is, therefore, diffeomorphic to the homogeneous space 
-4(2) x A(3)/H%. 

The families of curves Z^' C2 given by (8) have a large overlap. The following corollary elimi- 
nates this redundancy and, therefore, is useful for practical computations. 

Corollary 1 (reduced parallel projection criterion). A curve Z C M 3 projects to X C M? by 
a parallel projection if and only if there exist ai, a 2 , 6 £ R such that the curve X is -4(2)- equivalent 
to one of the following planar curves: 



z a u a 2 = {{z\ + aiz 3 ,z 2 + a 2 z 3 ) I (zi,z 2 ,z 3 ) £ Z), 



%b = {(zi + bz 2 ,z 3 ) | (zi,Z2,z 3 ) £ Z}, Z = {(z 2 ,z 3 ) \ (zi,z 2 ,z 3 ) £ Z). 



(10) 



Proof. We first prove that for any permutation (i, j, k) of numbers (1, 2, 3) such that i < j, and 
for any c±, c 2 £ R the set Z^ t ' C2 = {(z{ + c\Zk, zj + c 2 Zfc) | (z\, z 2; z 3 ) £ Z} is „4(2)-equivalent to 
one of the sets listed in (10). 
Obviously, Z,\^i 2 ■■ 

For Zcl 3 il, if c 2 ^ then 

equivalent to Z aiA2 with a\ = 
b = c\. 



ai>a2 with a\ = ci and a 2 = c 2 . 



3 i 
~ and a 2 



Zi + ClZ 2 
z 3 + c 2 z 2 



*1 - 1*3 

*2 + - 2 ZZ 



and so i?ci 3 c2 is -4(2)- 



-L. Otherwise, if c 2 = 0, the i?, 



1,3,2 
ei,C2 



^ with 
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Similarly for Zcl 3 cl, if C2 7^ then Zcucl is -4(2)-equivalent to Z ai;(l2 with a\ = ^ and 
02 = — Otherwise, if C2 = 0, then Zc^cl = [zi + c\Z\^Zz). If c% / then Zcf^ is A(2)- 

" 1 ~ 2 3 1 ~ 

equivalent to Z\, with b = —, otherwise c\ = and i? C i,c 2 = -2^ 

We can reverse the argument and show that any curve given by (10) is ,4(2)-equivalent to 
a curve from family (8). Then the reduced criteria follows from Theorem 2. ■ 

We note that the map 

x = z 1 +aiz 3 , y = Z2 + a 2 z 3 (11) 

is a parallel projection onto the z\, ^-coordinate plane in the direction of the vector (— ai, — 02, 1) 
with coordinates on the image plane induced from M 3 , namely, x = z\ and y = z%. We call (11) 
the canonical projection in the direction (— a%, —02, 1). The map x = z\ + bz2, y = z 3 is a pro- 
jection onto the z±, ^-coordinate plane in the direction of the vector (—6, 1,0) with coordinates 
on the image plane induced from M 3 , namely, x = Z\ and y = Z3, and finally the map x = z%, 
y = Z3 is the orthogonal projection onto the z%, Z3-plane. 



4 Solving the group-equivalence problem 

Theorems 1 and 2 reduce the projection problem to the problem of establishing group-action 
equivalence between a given curve and a curve from a certain family. In this section, we give 
a solution of the group-equivalence problem for planar algebraic curves. In Section 4.1, we con- 
sider a rational action of an arbitrary algebraic group on M 2 and define a notion of a classifying 
set of rational differential invariants. In Section 4.2, we define a notion of exceptional curves 
with respect to a classifying set of invariants and define signatures of non-exceptional curves. 
We then prove that signatures characterize the equivalence classes of non-exceptional curves. In 
Section 4.3, we produce explicit formulae for classifying sets of rational differential invariants for 
affine and projective groups. In Section 4.4, we specialize our signature construction to rational 
algebraic curves and provide examples of solving group-equivalence problem for such curves. 

We note that differential invariants have long been used for solving the group-equivalence 
problem for smooth curves. Classical differential invariants were obtained with the moving frame 
method [10], which most often produces non-rational invariants. Signatures based on classical 
differential invariants were introduced in [9]. For smooth curves, the equality of signatures of 
two curves implies that there are segments of two curves that are group-equivalent (in other 
words, these curves are locally equivalent), but the entire curves may be non-equivalent. This 
is well illustrated in [25]. In a recent work [21], a significantly more involved notion of the 
extended signature was introduced to solve global equivalence problem for smooth curves. 

The rigidity of irreducible algebraic curves allows us to use simpler signatures to establish 
global equivalence. Rationality of the invariants as well as explicit characterization of exceptional 
curves allows us to solve equivalence problem using standard computational algebra algorithms. 

4.1 Definition of a classifying set of rational differential invariants 

A rational action of an algebraic group G on IR 2 can be prolonged to an action on the n-th 
jet space J n = IR" +2 with coordinates (x,y,y^\ . . . ,y( n >) as follows 9 . For a fixed g G G, let 
(x, y) = g ■ (x,y). Then x, y are rational functions of (x, y) and 

g-{x,y,yW,...,y^) = (x,y,y^,...,y^), (12) 
9 Here y = y {0) and J° = R 2 . 
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where 



g W=#M -dfor .,„-! 5<* +1 > - S [S W (*,»,» (1) , • • • 



Here ^ is the total derivative, applied under assumption that y is function of x. 10 We note that 
a natural projection 7rj?: J n — > J k , k < n is equivariant with respect to action (12). For general 
theory of rational actions see [28] and for general definitions and properties of the jet bundle 
and prolongations of actions see [26]. 

Definition 5. A function on J n is called a differential function. The order of a differential 
function is the maximum value of k such that the function explicitly depends on the variable y^ . 
A differential function which is invariant under action (12) is called a differential invariant. 

Remark 4. Due to equivariant property of the projection tt2 : J n — > J k , k < n, a differential 
invariant of order k on J k can be viewed as a differential invariant on J n for all n> k. 

Definition 6. Let G act on M. N . We say that a set I of invariant rational functions on M. N 
separates orbits on a subset W C if W is contained in the domain of definition of each I £ I 
and Vwi,t«2 G W 

I(w\) = 1(1112) VI G 1 3g G G such that w\ = g ■ ui2- 

Definition 7 (classifying set of rational differential invariants). Let r-dimensional algebraic 
group G act on M 2 . Let K and T be rational differential invariants of orders r — 1 and r, 
respectively. The set X = {K, T} is called classifying if K separates orbits on a Zariski open 
subset W r ~ 1 C J 7 * -1 and X = {K, T} separates orbits on a Zariski open subset W r C J r . 

4.2 Jets of curves and signatures 

In this section, we assume that X C M? is an irreducible algebraic curve, different from a vertical 
line. Let F(x, y) be an irreducible polynomial, whose zero set equals to X . Then the derivatives 
of y with respect to x are rational functions on X, whose explicit formulae are obtained by 
implicit differentiation 

(1) _ _Fx (2) _ -F X X Fy + 2 F X y F X Fy - Fyy F% 

- F ' vx - F 3 , .... 

Definition 8. The n-th jet of a curve X C M 2 is a rational map jy : X > J n , where for p G X 

j x (p) = (x(p),y(p),y^(p),...,y^(p)). 

From the definition of the prolonged action (12), it follows that for all g G G and p G X the 
following equality holds, whenever both sides are defined. 

j g l .x(g-p) = g-[fx(p)]- (13) 

Definition 9. A restriction of a rational differential function J n — — > M to a curve X is 
a composition of with the n-th jet of curve, i.e. = $ o j^. If defined, such composition 
produces a rational function X > R. 



10 We note the duality of our view of variables y^ k \ On one hand, they are viewed as independent coordinate 
functions on J n . On the other hand, operator 4- is applied under assumption that y is a function of x and, 
therefore, j/W is also viewed as the k-th derivative of y with respect to x. (The same duality of view appears in 
calculus of variations.) 
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Definition 10. Let X = {K, T} be a classifying set of rational differential invariants for G-action 
(see Definition 7) . Then a point p G X is called X-regular if: 

(1) p is a non-singular point of X; 

(2) ^(pjer^audiKpjer; 

(3) dK 



/Oand 



J.X 



^0. 



An algebraic curve AT C M 2 is called non- exceptional with respect to I if all but a finite number 
of its points are X-regular. 

Lemma 1. Let I = {K,T} be a classifying set of rational differential invariants (see Defini- 
tion 7). Let X C R 2 be a non Z-exceptional curve defined by an irreducible implicit equation 
F(x, y) = 0. Then 

(1) K\x and T\x are rational functions on X and therefore there exist polynomials k\,k 2 G 
R[x, y] with no non-constant common factors modulo F, and polynomials ti,t 2 G R[x, y] 
with no non-constant common factors modulo F, such that 

K\x{x,y) = — r and T{x,y )\ x = — -. (14 

k 2 {x,y) t 2 {x,y) 

(2) The Zariski closure Sx of the image of the rational map S\ x : X ► R 2 , defined by 

S\x{p) = (K\x(jp)jT\x(p)) for p G X, is the variety of the elimination ideal X = X n 
R[x, t\, where 

X := (F, k 2 >c - ki,t 2 r - t±, k 2 t 2 o- - 1) C T,x,y, a}. (15) 

(3) dim Sx = if and only if Kx and Tx are constant functions on X and dim Sx = 1 
otherwise. In the latter case, Sx is an irreducible algebraic planar curve, i. e. a zero set of 
an irreducible polynomial Sx(k,t). 

Proof. (1) A function K\x '■ X >■ M is a composition of rational maps 

fx 1 : X- -> J r ~ x and K : J 1 " 1 - -> M 

(see Definitions 8 and 9). Since X is non-exceptional this composition is defined for all but 
finite number of points on X and therefore K\x is a rational function. The same argument 
shows that T\x is a rational function. Since X is defined by F(x, y) = 0, where F is irreducible, 
there exist polynomials k\,k 2 G M[x,y] with no non-constant common factors modulo F, and 
polynomials ti,t 2 G with no non-constant common factors modulo F, such that (14) 

holds. 

(2) By definition, 

{, ^ i -, / x x ki(x,y) ti(x,y) 

Sx = \ (k,t) I 3 (x,y) F{x,y) = A x = y 1 ^^ A r - 



k 2 (x,y) ' t 2 (x,y) 

and therefore is the projection of the variety defined by X, given by (15), to the x, r-plane. It is 
the standard theorem in the computational algebraic geometry (see, for instance, [12, Chapter 5]) 
that the Zariski closure Sx of this projection is the variety of the elimination ideal Xnl[x,r]. 

(3) It is not difficult to prove, in general, that the Zariski closure of the image of an irre- 
ducible variety under a rational map is an irreducible variety. The dimension of this closure is 
less or equal than the dimension of the original variety. Thus Sx is an irreducible variety of 
dimension zero when the signature map S\x = T\x) is a constant map and of dimension 

one otherwise. In the latter case Sx is an irreducible algebraic planar curve and, therefore, is 
a zero set of a single irreducible polynomial. ■ 
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Definition 11. Let X = {K, T} be a classifying set of rational differential invariants with 
respect to G-action and X be non-exceptional with respect to X. 

(1) The rational map S\ x : X > R 2 defined by S\x(p) = (K\ x (p),T\ x (p)) for p G X is 

called the signature map. 

(2) The image of S\x is called the signature of X and is denoted by S x . 

Theorem 3 (group-equivalence criterion). Assume that irreducible algebraic curves X\ and X<i 
are non- exceptional with respect to a classifying set of rational differential invariants I = (K, T) 
under the G-action. Then X\ and X2 are G-equivalent if and only if their signatures are equal: 

X\ = X 2 <;=>■ S Xl = S X2 . 
G 

Remark 5. Assume that two curves X\ and X2 have non-constant signature maps, and so the 
closures of their signatures are zero sets of polynomials S Xl (x, r) and S X2 (>c, r), respectively. 
The equality of signatures S Xl = S X2 , implies that S Xl (x, r) is equal up to a constant mul- 
tiple to S X2 (x,t). The converse is true over C, but not over K, because the latter is not an 
algebraically closed field (see Example 1 below and [12] for general results on implicitization). 

Proof of Theorem 3. Direction follows immediately from the definition of invariants. 
Below we prove <^=. We notice that there are two cases. Either K\ Xl and K\x 2 are constant 
maps on X\ and X2, respectively, and these maps take the same value. Otherwise both K\ Xl 
and K\ X2 are non-constant rational maps on X\ and X%, respectively. 

Case 1: There exists xq G K such that ^x^pi) = xq and Kx 2 (P2) = xq for all pi G X\ 
and for all p 2 £ X 2 . Since X\ and X2 are non-exceptional, we may fix I^-regular points 
Pi = (xi, y\) G X\ and p2 = (x2, 2/2) G X2. Then, due to separation property of the invariant K, 
3g G G such that j^T^pi) = g ■ [j^~ 1 (P2)]- We consider a new algebraic curve X% = g ■ X2. 
Then due to (13), we have 

^; 1 (pi)=^; 1 (pi)=:p ( ^ 1) - as) 

Since pi is a X-regular point of X%, it follows from (16) that it is also a X-regular point of X3 
and, in particular, is non-singular. Let F\(x, y) = and Fs(x, y) = be implicit equations of X\ 

pi 



and X3, respectively. We may assume that 7^ and 7^ (otherwise, 



and 9 -§ 



pi 



dy 



Pi 



7^ and we may use a similar argument). Then, there exist functions fi(x) and 

pi 

f^{x), analytic on an interval I 3 x%, such that F\(x, fi(x)) = and F^(x, fs(x)) = for x G I\. 
Functions y = fi(x) and y = fs(x) are local analytic solutions of differential equation 

K(x,y,yV,...,y^) = x (17) 

with the same initial condition /j \xi) = f^ k \xi), k = 0, ... ,r — 1 prescribed by (16). From 
the I- regularity of pi, we have that '" v 



yW=H(x,y,yW,...,y^), 



^ and so (17) can be solved for t/ r x ): 

p(r-i) r y ' y 



where function H is smooth in a neighborhood p( r_1 ) G J r_1 . From the uniqueness theorem for 
the solutions of ODEs, it follows that fi(x) = fs(x) on an interval I 3 x\. Since X\ and X3 are 
irreducible algebraic curves it follows that X\ = X3. Therefore, X\ = g ■ X%. 

Case 2: K\ Xl and K \ X2 are non-constant rational maps. Then S Xl = S X2 is a one-dimensional 
set that we will denote S. Let S(x, r) = be the implicit equation for S (see Lemma 1). We know 
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that ff (x, t) 7^ for all but finite number of values (x, r), because, otherwise, Jf| < v t and K\x 2 
are constant maps. Therefore, since the curves are non-exceptional, there exists X-regular points 
Pi = (%i,yi) G X% and p2 = (x2,U2) G X% such that 



^I-Vi(pi) = K\x 2 (p 2 ) =: x , ^(^(pi) = T\x 2 (p 2 ) =: r , 



as 



(x ,r )/0. (18) 



Due to separation property of the set Xq = {K,T}, 3g £ G such that j x (pi) = g • [jx 2 (P2) 
We consider a new algebraic curve X% = g ■ X 2 . Then due to (13), we have 



From (18), (19) and Z-regularity of the point pi G X\ it follows that 

dT 



tf(pM) 



r( P W) 



and 



<9yM 



,(r) 



^0. 



(19) 



(20) 



Since pi is a X-regular point of A4, it follows from (19) that it is also a X-regular point of X$ 
and, in particular, is non-singular. Let F\(x, y) = and F^(x, y) = be implicit equations of A4 



and X3, respectively. We may assume that 7^ and ^j- 7^ (otherwise, 



and^ 



pi 



dy 



Pi 



7^0 



pi 



Pi 



7^ and we may use a similar argument). Then, there exist functions fi(x) and fs{x), 

analytic on an interval / 3 xi, such that F\(x, fi(x)) = and F^(x, fs(x)) = for x G Ji. Then 
functions y = /i(cc) and y = /s(x) are local analytic solutions of differential equation 



S(K{x, y, y « . . . , y^ 1 )) , T(x, y, y« . . . , y^)) 







(21) 



with the same initial condition f[ k \xi) = f^'(xi), k = 0, . . . ,r, dictated by (19). 



ST 



Since w (x ,t ) + and 
for y( r ^: 



,(«0 



7^ (see (18) and (20)), equation (21) can be solved 



,{r) 



H(x,y,y 



(i) 



,(r-l) 



)• 



where function is smooth in a neighborhood p( r ) G J r . From the uniqueness theorem for 
the solutions of ODE it follows that fi(x) = fo{x) on an interval I 3 x\. Since X\ and X3 are 
irreducible algebraic curves it follows that X\ = X3. Therefore, X\ = g ■ X 2 . ■ 

From the proof of Theorem 3 we may deduce the following: 

Corollary 2. 

1) Assume K\x 1 is a constant function with value xq G R, i.e. X|^ 1 (p) = >cq for all p G X\. 
Then S\x 2 = S\x 1 if and only if K\x 2 (p) = xo f or M P G X 2 . 

2) //dim(«S|^ nS|* 2 ) = 1 then S\ Xl = S\x 2 - 

4.3 Classifying sets of invariants for affine and projective actions 

In this section, we construct a classifying set of rational differential invariants for affine and 
projective actions. We will build them from classical invariants from differential geometry [4, 10]. 
We start with Euclidean curvature 



y 



(2) 



;i + [y (l)]2)3/2' 



(22) 
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which is, up to a sign 11 , a Euclidean differential invariant of the lowest order. Higher order 
Euclidean differential invariants are obtained by differentiating the curvature with respect to 
the Euclidean arclength ds = y/l + [yW] 2 dx, i.e. k s = «g = -j=L= ^, k„ = 

Equi-affine and projective curvatures and infinitesimal arclengths are well known, and can 
be expressed in terms of Euclidean invariants [13, 24]. In particular, 54.(2)-curvature p and 
infinitesimal 54.(2)-arclength da are expressed in terms of their Euclidean counterparts as fol- 
lows 

3k(k ss + 3k 3 ) — 5k 2 i/o 
p = -7- , da = k ' as. 

I" 9k 8/3 

By considering effects of scalings and reflections on 54.(2)-invariants, we obtain two lowest order 
A ( 2 ) - invar i ant s 



K A = ^f, T4=^f. (23) 



0«) 2 T A = - 

p 3 ' 3/i 

They are of order 5 and 6, respectively, and are rational functions in jet variables. 

■p<y£(3)-curvature r] and infinitesimal arclength dp are expressed in terms of their SA- 
counterparts 

_ — 7p aa — 9 fig /J, . _ 1/3 , 

'/ — g ; 3 j U P — Pa 

6p a ' 

The two lowest order rational 'P^£(3)-invariants are of differential order 7 and 8, respectively 

K v = r/ 3 , T v = n p . (24) 

Explicit formulae for invariants in terms of jet coordinate are given by (46) and (47), (48) in 
Appendix A. 

Theorem 4. According to Definition 7: 

(1) The setX A = {K A ,T A } given by (23) is classifying for the A(2)-action on I 2 . 

(2) The setT-pgc = {Kp,Tp} given by (24) is classifying for the VQC(3) -action onM?. 
Proof. We start by introducing differential functions 

A i: =3 y ( 4 ) y ( 2 )-5[ y ( 3 )] 2 , (25) 
A 2 :=9y( 5 )[y( 2 )] 2 -45y( 4 )y( 3 V 2 )+40[y( 3 )] 3 . (26) 

(1) We note that dim .4(2) = 6. We will prove the separation property of Z4 = {K A ,T A } 
given by (46) on a Zariski open subset 

W 6 = {p {6) E J 6 I y {2) and Ai + 0} 

of J 6 , where Ai is given by (25), and the separation property of K A on W 5 = tt®(W 6 ) C J 5 . 
An affine transformation can be written as a product of a Euclidean, an upper triangular area 
preserving linear transformation, a scaling and a reflection 



.4 




/ 
1 




where c 2 + s 2 = 1, e = ±1 and h / 0. 

11 The sign of k changes when a curve is reflected, rotated by ir radians or traced in the opposite direction. 
A rational function k 2 is invariant under the full Euclidean group. 
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A Euclidean transformation with 

1 t/ 1 ) yW y + x t/ 1 ) x — y 



Vi + b (1) ] 2 ' ' \A + b (1) ] 2 ' Vi + b (1) ] 2 ' Vi + b (1) ] 2 ' 

brings a point = (x,y,y^ ,y^ ,y^ ,y( 4 \y^) to p^ = (0, 0, 0, y^\ y^\ y[^\ y^) , where 
y^ if and only if y( 2 ) 0. An upper triangular transformation with e = [yj 2 ^] l ^ and 

/ = 3 [^2)|5/3 brings the point to a point p^ = (0, 0, 0, 1, 0, y^\y^) , where ^ if and 
only if 3y( 2 )y( 4 ) — 5[y( 3 )] 2 ^ 0. A scaling transformation with h = [y^] 5 brings to a point 
Pg = (0, 0, 0, 1, 0, 1, y 3 ) ■ Finally, if y 3 < 0, we can apply a reflection (e = —1) to bring p^ 
to pi 5) = (0, 0, 0, 1,0, l,2/i 5) ) with yf ] > 0. 

By the above argument, for any p( 5 ),q( 5 ) E W 5 , there exist A, B E A(2) such that p[ 5 ^ = 
A-p( 5 ) = (0,0,0, 1,0, l,yf } ) and qf } = B ■ q( 5 ) = (0, 0, 0, 1, 0, 1, y^), where yf } > and 
2/2 — 0- Then, since -fT^ is invariant, 

^(P (5) ) = K A [jg>) = 3[y{ 5) ] 2 and K A (^) = K A [<g>) = 3[yf] 2 . 

Since y[ 5 ^ > and y^ > 

^( P (5) )=^(q (5) ) ^ yf , = yi ) <=► P S 5) = qi 5) ^ q ( 5 ) = i?A- 1 .p( 5 ), 

and, therefore, K A is separating on W 5 . 

Similarly, for any p (6) ,q (6) E W e , there exist A, B E .4(2) such that pp = vl • p (5) = 
(0,0,0,l,0,l,y{ 5) ,yi 6) ) and qf } = B ■ = (0, 0, 0, 1, 0, 1, yf , yf ) , where yf } > and 
2/2 > 0. Then since and T4 are invariant 

^(P (6) ) = ^(PS 6) ) = 3[y$ 5) ] 2 , T^(p( 6 )) = T4(pf)) = - 5, 

and 

^(q (6) ) = K A {$) = 3[yf] 2 , T^(q( 6 )) = 3U(qW) = yf " 5- 

Then 

Mp W )=Mq W ), T^(p( 6 ))=^(q( 6 )) o yf)=yf, y« = yf 
O pi 6) =qS 6) «> q( 8 )=Bi4- 1 .pW 

and thus the separation property of {if^T^} on J 6 is proved. 

(2) We note that dimVQC(3) = 8. We will prove the separation property of X-p = {Kp,Tp} 
given by (47) and (48) on a Zariski open subset 

W 8 = {p {8) E J 8 I y (2) ^ and A 2 / 0} 

of J 8 , where A2 is given by (45), and the separation property of K-p on W = 7r 8 (VF 8 ) E J 7 . 
Observe that any element in A E VQC(3) can be written as the following product 



.4 




/ 
1 

1 
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where c 2 + s 2 = 1, e / and g ^ 0. In the first part of the proof, we have shown that the two 
matrices on the right can bring a point p( 8 ) G W 8 to the point = (0, 0, 0, 1, 0, y[^\ . . . , y^) ■ 
Observe that W 8 is invariant under these transformations. Therefore p^ G W 8 and hence 

(5) 

y\ 7^ 0. The last transformation with 



a = m 

and 



(5)1 1/3 



5[^ 4) ] 2 - yj 6) 



._ [yf ] ] 2 -iQyf ) bi 4) ] 2 -3bi 5) ]M 4) + 25[^] 4 

18^] 2 

brings p^ to = (0, 0,0,1, 0, 0, 1, 0, y£ , V^) ■ A direct computation shows that 
Mp? } ) = bfT and T v (p { ?)=3yi 8) - 1 f. 

(7) (S) 

We observe that y^ and y\ can be uniquely determined from the values of invariants K-p 
and T-p and therefore we can complete the proof of the separation property of Kp on W 7 and 
the separation property of {Kp,Tp} on W 8 by an argument similar to the one presented in 
part 1 of the proof. ■ 

Theorem 4, in combination with Theorem 3, leads to a solution for the projective and the 
affine equivalence problems for non-exceptional curves. The following proposition describes 
exceptional curves. 

Proposition 2. X4- and Xpgc- exceptional curves are lines and conies 12 . 

Proof. In the affine case, we note that j\{p) G J 5 \ W 5 and j%{p) G J 6 \ W 6 if an only if 
k Iat(p) = or fi\x(p) = 0. If k|^(p) = for more than finite number of points on X, then X 
is a line, if (x\x(p) = for more than finite number of points on X then it is a parabola (see 
Proposition 3). From the explicit formulae (46) we see that 



8K A 



dy&) 



18A 2 [ y ( 2 )] 2 8T A 



,(■)" A? ^ By® 



9[y 



(2)1 



,(«) 



A? 



These rational functions are not identically zero on X if neither = nor y^\x = 0, or 

equivalently X is not a line or a conic. Therefore, if X is not a line or a conic, it is {/C^,7yi}- 
regular. In the projective case, we note that jx(p) ^ J 7 \ W 7 and j%(p) G J 8 \ W s if an only if 
k |a , (p) = or fi, a \x(p) = 0. If fj, a \x(p) = for more than finite number of points on X then X 
is a conic (see Proposition 3). From the explicit formulae (47) and (48), we see that 



dK v 



dyW 



, ( 7) dyW 



7^0 

,(8) 



for all p( 7 ) G W 7 and all p( 8 ) G W 8 . Therefore, if an algebraic curve is not a line or a conic it 
is {/Cp, 7p}-regular. ■ 



12 By a conic we mean an irreducible real algebraic planar curve of degree 2, i.e. a parabola, an ellipse, or 
a hyperbola. 
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Remark 6. It is well known (and easy to prove) that the set of all lines constitutes a single 
equivalence class (an orbit) under both A(2, M)-action and VQC{3, R)-action. Under the A(2, In- 
action the set of all conies splits into three orbits: the set of all parabolas, the set of all hyperbolas 
and the set of all ellipses, although under the A(2, C)-action the set of all hyperbolas and ellipses 
comprise a single orbit. All conies constitute a single orbit under VQC{3, M)-action, see [3, 
Section II. 5]. Therefore an X_4-exceptional algebraic curve is not A(2, C)-equivalent to a non Z4- 
exceptional algebraic curve and an 2/pg£-exceptional algebraic curve is not VQC(3, C)-equivalent 
to a non X-pg^-exceptional algebraic curve. 

The projective and the affine equivalence problems for exceptional curves can be easily solved 
using the above remark and the following proposition. 

Proposition 3. Let X be an irreducible planar algebraic curve. Then 

(1) X is a line <^=^ k\x = 0; 

(2) X is a parabola ^=>- \i\x = 0; 

(3) X is a conic <J=^ fi a \x = 0, 

where "= 0" means that a corresponding rational function is zero at every point of X . The above 
statements are true for both real and complex algebraic curves. If X is a real algebraic curve, 
then it is a hyperbola if and only if y\x is a negative constant, while X is an ellipse if and only 
if \i\x is a positive constant. 

The proof of part (1) of Proposition 3 follows immediately from (22). Proofs of the other 
statements can be found in [18, Section 7.3]. The following corollary is obtained from Proposi- 
tion 3 using explicit formulae for equi-affine invariants 

_ 1 Ai 1 A 2 

M ~9[ y (3)]8/s' ^-27[ y (2)f 

where Ai and A2 are given by (25) and (26). 

Corollary 3. Let X be an irreducible planar algebraic curve. Assume that X is not a vertical 
line. Let Ai and A2 be given by (25) and (26). Then 

(1) The restrictions A±\x and A2\x ar & rational functions on X. 

(2) Ai\x is a zero function if and only if X is a line or a parabola. Otherwise, restrictions of 
A(2) -invariants Kj\x andTj\x are rational functions of X , 13 

(3) A%\x is a zero function if and only if X is a line or a conic. Otherwise K-p\x and Tp\x 
are rational functions on X. 



4.4 Signatures of rational curves 

In this section, we adapt the signature constructions to rational algebraic curves and give 
examples of solving the affine and projective equivalence problems using signatures. We adapt 
Definition 9 to rational curves as follows. Let X is a rational curve parameterized by j(t) = 
(x(t),y(t)), such that x(t) is not a constant function 14 . Make a recursive definition of the 
following rational functions of t: 

X X 



13 Ka\x and Ta\x are defined and are both zero functions when X is either an ellipse or a hyperbola. We know, 
however, that ellipses and hyperbolas are not A(2, R)-equivalent. There is no contradiction with Theorem 3, 
because ellipses and hyperbolas are Z^-exceptional per Definition 10. 

14 Equivalently, X is not a vertical line. 
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Figure 2. X\\ x A + if - lOxy = 0, X%: y - xy + 1 = 0, A' 3 : y A - x 2 + xy 2 = 0. 



where ' denotes the derivative with respect to the parameter. Let $ be a rational differential 
function. Then the restriction of 3>| 7 is computed by substituting (27) into <3>. If defined, $L is 
a rational function of t. Recalling Definition 11 of signature and Corollary 3, we conclude that: 

Proposition 4. Let X be an irreducible planar algebraic curve parameterized by a rational 
map 7(t). Assume A2| 7 7^ ; where A2 is given by (26). Then 



(1) the A(2) -signature of X is the image of a rational map 5^| 7 (t) := (-?Gl| 7 , TU| 7 ) from 



to 



(2) the VQC(3)- signature of X is the image of a rational map S*p| 7 (t) := (X-p| 7 , Tp| 7 ) /) 



rom . 



to 



The signatures can be computed either using inductive formulae for invariants given by (23) 
and (24), or explicit formulae given by (46) and (47), (48) in Appendix A. A solution for A(2)- 
and "P<5£(3)-equivalence problems for rational curves follows from Theorem 3 and is illustrated 
by the following examples. Invariants and signatures were computed using Maple code [7]. 

Example 1 ("P^£(3)-equivalence problems). Consider three rational cubics pictured on Fig. 2. 
These curves have the following rational parameterizations: 



X\ is parameterized by a(t) = 
X2 is parameterized by /3(s) = 
X3 is parameterized by j(w) 



Wt Wt 2 



t 3 + 1 ' t 3 + 1 
s 3 + 3s 2 + 3s + 2 



w" 



s + l 



u>- 



-,s + l 



v u> + 1 w + 1 / 

The projective signature of X\ is parameterized by invariants 



K v \ a {t) 



9261 (t 6 - t 3 + 1)H 3 
^0 (t 3 - l) 8 : 



TvUit) 



21 (t 3 + l) 4 
10 (i 3 - 1) 4: 



(28) 
(29) 
(30) 

(31) 



while the signature of X% is parameterized by invariants 

9261 (s + l) 3 (s 6 + 6s 5 + 15s 4 + 19s 3 + 12s 2 + 3s + l) 3 



KvHs) 



50 (s 2 + 3s + 3) 8 s 8 

21 (s 3 + 3s 2 + 3s + 2) 4 
10 (s 2 + 3s + 3) 4 s 4 



(32) 
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Figure 3. X 4 : y 3 



xy 2 + x 2 



xy + 2x-3y + 2 = Q. 



Although it is not obvious, the curves defined by parameterizations (31) and (32) satisfy the 
same implicit equation 



62523502209 + 39697461720r - 6401203200x + 5250987000r 2 - 2032128000xr 

0. (33) 



+ 163840000x 2 + 259308000r 3 + 53760000xr 2 + 4410000t 4 



This is a sufficient condition for the equality of signatures Sx x and Sx 2 over complex numbers, 
but it is not sufficient over reals. We can look for a real rational reparameterization t = 4>(s) 
by solving a system of two equations K-p\ a (t) = K-p\p(s) and T-p\ a {t) = T-p\p(s) for t in terms 
of s. One can check that t = s + 1 provides a desired reparameterization. Thus Sx x = Sx 2 and 
hence, by Theorem 3 

Xi = X 2 . 

Reparameterization t = s + 1 allows us to find pairs of points on X\ and X 2 which can be 
transformed to each other by Pt/£(3)-transformation that brings X\ to X 2 - Since four of such 
pairs in generic position uniquely determines a transformation we can compute that X 2 can be 
transformed to X\ by a transformation 



lOy 

■ 

X 



y 



10 

x 



It turns out that cubic X% has constant 'P<5£(3)-invariants 
250047 



K p \ 7 (w) 



12800 



and T-p\j(w) = 



(34) 



(35) 



and therefore its signature degenerates to a point. Thus, by Theorem 3, X3 is not VQC{?>)- 
equivalent to either X\ or X 2 . 

To underscore the difference between the solution of the equivalence problems over real and 
over complex numbers, we will consider one more cubic X4 pictured on Fig. 3, whose rational 
parameterization is given by 



S(u) 



u 6 + 1 u 2 + 1 



II 



II 



The signature of X4 is parameterized by invariants 



K v \ s (u) 



6751269 (u - l) 3 (u 2 -u + lf(u 3 



3u 2 + 1) V 



50 



it" 



6u 2 + 3u + 1) 
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189 (u 3 + 1 - 3u) 4 

T v \s{u) - 



10 (n 3 - 6n 2 + 3n + l) 4 ' 

Invariants K-p\g and Tp\g satisfy the implicit equation (33). Since the signatures of X\ and X 2 
satisfy the same implicit equation, we can conclude that, over the complex numbers, X4 is 
projectively equivalent to both X\ and X2. In fact, we can find that for 

1 1 ^ 3 
V3i 

2 2 

(where we are free to choose any of the three cubic roots) the complex projective transformation 

-6iV3ry + 3r 2 (iy/3 - 3) _ _ 2x - r (5iV3 + S)y + r 2 (iy/3 - 9) 

X ~ 2x - 2iy/3ry + r 2 (iy/3 - 3) ' 2x - 2iy/3ry + r 2 (iy/3 - 3) 

transforms X4 to X 2 . Our attempt to solve Tp\$(u) = Tp\p(s) and Tp\g(u) = Tp\p(s) for u in 
terms of s gives a rather involved rational complex reparameterization that transforms signature 
map of 5 into the signature map of /?, but no real reparameterization. Therefore 

X 4 ^ X 2 but X 4 % X 2 . 

vgc{3,c) vgc{3,m) 

Example 2 („4(2)-equivalence problems). We can again consider three cubics pictured on Fig. 2, 
but now ask if they are „4(2)-equivalent. Recalling that .4(2) is a subgroup of VGC(S) we can 
immediately conclude from the previous example that X3 is not 4.(2)-equivalent to either X\ 
or X 2 . To resolve the equivalence problem for X\ and X 2 we need to compute their affine 
signatures. The affine signatures of X\ is parameterized by invariants 

1 (t 3 - I) 2 1 
KMt) = -^^T^' = -Jo' 

It turns out, that restrictions of both invariants, Kj\\p{s) and Tj\\p(s), are non-constant functions 
of s. Hence, X\ and X 2 have different affine signatures. Therefore 

X\ ¥ X 2 . 

A{2) 

In fact, affine signatures for all four curves X\, X 2 , X3 and X^ have different implicit equations 
and therefore no two of them are affine equivalent neither over real numbers, nor over complex 
numbers. 



5 Algorithm and examples 

The algorithms for solving projection problems are based on a combination of the projection 
criteria of Section 3 and the group equivalence criterion of Section 4. 

5.1 Central projections 

The following algorithm is based on the central projection criterion stated in Theorem 1 and 
the group-equivalence criterion stated in Theorem 3. In the algorithm, we compute restrictions 
of differential functions A2, K-p and T-p to a curve parameterized by j(t) and to a family of 
curves parameterized by e(c, s), where c = (ci, c 2 , C3) determines a member of the family and s 
serves to parameterize a curve in the family. These restrictions are computed by substitution 
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of (27) into formula (26) for A2 and into (47) and (48) for K-p and Tp, respectively. When the 
restrictions to e(c, s) are computed, derivatives in (27) are taken with respect to s. 

One can use general real quantifier elimination packages, such as Reduce package in Mathe- 
MATICA, to perform the steps involving real quantifier elimination problems. To make an 
efficient implementation, one needs to take into account specifics of the problems at hand. 
This lies outside of the scope of the current paper and is a subject of our future work. 

Algorithm 1 (central projections). 

Input: Parameterizations F = (z±, z%, Z3) G Q(s) 3 and 7 = (x, y) G Q(t) 2 of rational algebraic 
curves Zcl 3 and X C M. 2 , respectively, such that Z is not a line, that is, T x T / 0. 

Output: The truth of the statement: 



3 [P] G CV, such that X = P{Z). 
Steps: 

1. [If X is a line, then determine whether Z is coplanar.] 



if 



m) 



0, then return the truth of the statement 



») 



0. 



2. [Describe a family of parametric curves where c = (c\, C2, C3) specifies a member of the family.] 

-= (S^)eQ( Cl ,C 2 ,C3, S ) 2 . 

3. [Compute the denominators of the rational invariants.] 

compute A2| 7 and A 2 | e , using the formula (45) in Appendix A. 

4. [If X is a conic, then determine whether, for some c, the rational map e(c, s) parameterizes 
a conic] 

if A2I7 = 0, then return the truth of the statement 



3cG 



Z3 + c 3 / OA 

R(s) 



^ A A 2 | 

R(s) 



(•) 



0. 



(36) 



5. [If Z cannot be projected to a curve of degree greater than 2, then return FALSE.] 
if A 2 | e = 0, then return false. 

R(c,s) 

6. [Compute the rational invariants.] 

compute K-p\-y, Kp\ e , T-p| 7 and T-p\ e using the formulae (47), (48) in Appendix A. 

7. [Determine whether, for some c, the signature of the Zariski closure Z c of the curve parame 
terized by e(c, s) equals to the signature of X .] 

if K-p\-y is a constant rational function, then return the truth of the statement 
3 c G M 3 



c is generic A K-p \ e = K-p \ 

R(s) 



7' 



(37) 



else return the truth of the statement 

3 c G M 3 : c is generic A Kp\ e is not a constant rational function, 
AVsGM A 2 | e /0 3 i G M Kp\ e = Kp\ y A Tp\ e = Tp\ 7 , 



(38) 



where we define 
[c is generic] 



Z3+C3 + OA 
R(s) 



/OA A 2 | £ / 

R(s) R(s) 



(39) 
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Proof of Algorithm 1. On the first step of Algorithm 1, we consider the case when X is 
a line. Then Z can be projected to X if and only if Z is coplanar. Both conditions can be 
checked by computing determinants of certain matrices. If X is not a line we define, on Step 2, 
a rational map e that parameterizes a family of curves. On Step 3, we compute restrictions 
of differential function A2 to j(t) and e(c, s). We remind the reader, that in the latter case 
derivatives are taken with respect to s. Since on this step we know that X is not a line and 
there are values of c for which Z c is not a line, these restrictions are defined. On Step 4, we 

consider the case when X is a conic. Equivalently, by Corollary 3, A2U = 0. Then Z can be 

Q(t) 

projected to X if and only if 3 c such that the Zariski closure Z c of the curve parameterized 
by e(c, s) is projectively equivalent to X and therefore is a conic (see Remark 6). Equivalently, 
A2| e = 0. If X is not a conic, we reach Step 5, where we check for a possibility that Z c is 

M(s) 

a conic for all parameters c (equivalently A2| e = 0) and therefore it can not be projected 

K(c,s) 

to X, which at this step is known to have a higher degree. If that is not the case, we proceed to 
Step 6, where we compute restrictions of differential invariants K-p and T-p to 7(f) and e(c,s). 
Since on this step we know that X is of degree greater than 2 and there are values of c for 
which Z c is of degree greater than 2, these restrictions are defined. 

On Step 7, where we know that X is non-exceptional and decide if there exists c G M 3 
such that: (1) Z c is non-exceptional which is equivalent to condition (39); (2) the signatures 
of the algebraic curves X and Z c are the same. On Step 6, we computed rational functions 
K-p\ e (c, s) and Tp\ t (c, s). We need to show that if we substitute a specific value c = Co G M 3 
into these functions we obtain the same rational functions of s as we would obtain by com- 
puting restrictions of the invariants to the curve parameterized by e(co,s). For generic values 
of c defined by (39), we can show that this is true. Indeed, it is well known that taking 
derivatives with respect to one of the variables and specialization of other variables are com- 
mutative operations. From condition (39) it follows that the curve Z c is not a line and so 
the denominators of (27) are not annihilated by such specialization. Therefore the restric- 
tion of jet variables to a curve parameterized by e(c, s) commutes with a specialization of c. 
From condition (39) it also follows that the denominators of (47) and (48) are not annihi- 
lated by a generic specialization. Therefore, for a generic Co, rational functions Kp\ t (co,s) and 
Tp\ e (co,s) equal to the restrictions of the invariants K-p and T-p to e(co,s). To decide equality 
of signatures of the algebraic curve parameterized by 7(t) and the curve Z c we use Corollary 2 
with (37) analyzing the case of constant invariant K-p\x and (38) the case of non-constant 
invariant Kp\x- I 

Remark 7 (reconstruction). If the output is true, then, in many cases, in addition to es- 
tablishing the existence of c\, c%, C3 in Step 4 or 7 of the Algorithm 1, we can find at least 
one of such triplets explicitly. We then know that Z can be projected to X by a projection 
centered at (— c\, — C2, — C3). We can also, in many cases, determine explicitly a transformation 
LA] G VQC{2>) that maps X to the Zariski closure Z c of the image of the map e(c, s). We then 
know that Z can be projected to X by the projection [P] = [A] [P®] [B] , where P® and B are 
defined by (5). 

Example 3. We would like to decide if the spatial curve Z, pictured on Fig. 4 and parameteri- 
zed by 

r(s) = (s 3 ,s 2 ,s), seR (40) 

projects to any of the four cubics planar cubics described in Example 1. 

We start with cubics X±, X2 pictured on Fig. 2, whose parameterizations are given by (28) 
and (29), respectively. Since these two cubics are PC/£(3)-equivalent then the twisted cubic can 
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be either projected to both of them or to none of them (see Proposition 1). Let us "run" the 
algorithm for X2. Following Algorithm 1, since X^ is not a line, we proceed to Step 2 and define 
a family of curves 

e(ci,c 2 ,C3,s) = ( S +Cl , S ^ C2 V (41) 
V s + c 3 s + c 3 J 

Since Xi is not a conic, we proceed to Step 6. Invariants for X\ and Xi are given by formulae (31) 
and (32), respectively. The formulae for rational functions K-p\ e {c\,C2, C3, s) and Tp\ e (ci,C2, c 3 , s) 
are too long to be included in the paper, but can be computed using our Maple code [7]. Step 7 
of Algorithm 1 returns true. In fact, c = (1,0,0) satisfies conditions (38) and therefore Z can 
be projected to both X\ and X% with a camera centered at (—1, 0, 0). We note that the canonical 
projection 

x = ^, y= 23 . (42) 

Z2 Z 2 

with this center maps Z to X2 (although resulting parameterization differs from (29).) From 
decomposition (4) we know that a projection from Z to X\ is the composition of the projec- 
tion (42) and VGC(3) transformation (34) that brings X2 to X\. The resulting projection is 

rr - 10£3_ „ _ 10^2 

x — 21+1 >y — zi+i- 

As a side remark, we note K-p\ e (c, s) and Tp\ e (c, s) satisfy the implicit equation(33) indepen- 
dently of c. Thus, the closure of the signature Z c belongs to the closure of the signature of X2 
(and X\). Therefore, we could immediately establish existence of a complex projection, because 

the conditions (39) for c being non-exceptional and such that K-p\ e 7^ define a Zariski open 

C(«) 

non-empty subset of C 3 . 

Let us now consider X^ pictured on Fig. 3. Again Algorithm 1 reaches Step 7 and returns 
true. In fact, c = (1, 1, 0) satisfies conditions (38) and therefore Z can be projected to X^ with 
a camera centered at (—1, —1, 0). Indeed, the canonical projection with this center 

21 + 1 Z 2 + \ 

x = , y = 

Z3 z 3 

maps Z to X4. 

This is a good point to compare real and complex projection problems. As was established in 
Example 1, the signatures of all three curves X\, X2 and X4 have the same implicit equation (33). 
Therefore, they are all VQC(3, C) -equivalent. Therefore, since there is a projection from Z to X\ 



24 



J.M. Burdis, LA. Kogan and H. Hong 



and X2, centered at (—1, 0, 0), there is a complex projection centered at (—1, 0, 0) from Z to X4. 
We also established, in Example 1 that 

X 4 ¥ X 2 (43) 
P6£(3,R) 

and, therefore, there is no real projection centered at (—1,0,0) from Z to X^ which, as we have 
seen, does not preclude the existence of a real projection with a different center (—1, —1, 0). 

Finally we consider X$, pictured on Fig. 2 with parameterization (30). From Example 1 we 
know that invariants for X3 are constants, see (35). Following Algorithm 1, we need to decide 
whether there exists c 6 M 3 , such that e(c, s) does not parameterize a line or a conic and 

K vU o,s)= 250047 V S £». 
' 1 V ; 12800 

This is, indeed, true for a = C2 = and C3 = 1. This is sufficient to conclude the existence 
of a real projection. We can check that Z can be projected to X3 by the a central projection 



23 + I ' » z 3 + l- 

The above example underscores Remark 1: although the twisted cubic can be projected to 
each of the planar curves X±, X2, X3 and X$, the planar curve X3 is not VGC(3, C)-equivalent 
to Xi, or X2, or A4. Also X4 is not P^£(3, IR)-equivalent to X\ or A^. 

Example 4. In this example, we establish that the twisted cubic (40) can be projected to any 
conic. As before, the family of rational curves e is defined by (41). The twisted cubic can be 
projected to a conic if and only if there exists c such that e(c, s) does not parameterize a line and 

A.2 1 e = 0. We can easily check that e(c, s) is not a line for all c and that A2\ t = whenever 

C(s) CO) 

ci = 4, c 2 = -cl (44) 

Let C3 = —a, where a is any real number. From (44) it follows that e parameterizes a conic 
if and only if c = (—a 3 , —a 2 , —a). The corresponding canonical projection centered (a 3 , a 2 , a) 
maps the twisted cubic to the parabola 

y + ay — x + a =0. 

Since all conies are 'P^£(3)-equivalent, we established that the twisted cubic can be projected 
to any conic. Moreover, we established that the twisted cubic is projected to a conic if and only 
if the center of the projection lies on the twisted cubic. 

So far in all our examples the outcome of the projection algorithm was true. Below is an 
example with false outcome. 

Example 5. We will show that the twisted cubic (40) can not be projected to the quintic 
co(t) = (t,t 5 ). The signature of the quintic is parameterized by a constant map 

1029 

Kr\ u {t) = and Tp\ u (t) = 0, Vt 

Following Algorithm 1, we need to decide whether there exists c£K 3 , such that e(c, s) does not 
parameterize a line or a conic and 

1029 

K v \ e (c,s) = —, VseE. 

Substitution of several values of s in the above equation yields a system of polynomial equations 
for 01,02,03 G K that has no solutions. We conclude that there is no central projection from Z 
to uj(t) = (i,i 5 ). This outcome is, of course, expected, because a cubic can not be projected to 
a curve of degree higher than 3. 
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5.2 Parallel projections 

The algorithm for parallel projections is based on the reduced parallel projection criterion stated 
in Corollary 1. This algorithm follows the same logic but has more steps than Algorithm 1, be- 
cause we need to decide whether a given planar curve is .A(2)-equivalent to a curve parameterized 
by a(s) = (z2(s), z 3 (s)), or to a curve parameterized by (3(b, s) = (zi(s) + bz2(s), z 3 (s)) for some 
b € M, or to a curve parameterized by 5(ai,a,2,s) = (z\(s) + 01^3(5), Z2 + «2 ^(s)) for some 
a = (01,02) G M 2 . Since the affine transformations are considered, projective invariants are 
replaced with affine invariants (see (46)). Due to its similarity to Algorithm 1, we refrain from 
writing out the steps of the parallel projection algorithm and content ourselves with presenting 
examples. A Maple implementation of the parallel projection algorithm over complex numbers 
is included in [7]. 

Example 6. As a follow-up to Example 3, it is natural to ask whether the twisted cubic can be 
projected to any of the cubics considered in that example by a parallel projection. Our imple- 
mentation of the parallel projection algorithms [7] provides a negative answer to this question, 
the twisted cubic can not be projected to X\, or X2 or X3, or X4 under a parallel projection 
even over complex numbers. There are plenty of rational cubics to which the twisted cubic can 
be projected by a parallel projection. For example, the orthogonal projection to the z\, Z3-plane 
projects the twisted cubic to (s 3 , s). 

Example 7. As a follow-up to Example 4, we consider the problem of the parallel projection 
of the twisted cubic to a conic. To answer this question we first consider a curve a(s) = 
(z2(s), z 3 (s)) = (s 2 , s), which is a parabola, and therefore the twisted cubic can be projected to 
any parabola. We then define a one-parametric family of curves 

f3(b, s) = ( Zl (s) + bz 2 (s), z 3 (s)) = (s 3 + bs 2 , s) 

and a two-parametric family 

8(a, s) = (zi(s) + aiz 3 (s), z 2 + a 2 zs(s)) = (s 3 + ais, s 2 + a 2 s). 

Obviously, there are no parameters o or b such that a curve in those two families becomes a conic 
(we can check this formally by computing rational functions A2\ a and A2I/3 and seeing that there 
are no values of a or b that will make them zero functions of s). Thus we conclude that the 
twisted cubic can not be projected to either a hyperbola or an ellipse by a parallel projection. 

For a less obvious example and to finally get away from the twisted cubic we consider the 
following: 

Example 8. We would like to decide whether the spatial curve Z parameterized by 

r(a) = (s 4 ,s 2 ,s), seR 
can be projected to X parameterized by 
j(t) = (t 4 + t 3 ,t 2 + t), tel. 
The signature of X is parameterized by invariants 

(24i 5 + 51f 4 + 57t 3 + 33t 2 + 9t + l) 2 



K A \ 7 (t) = -1600 
T4| 7 (t) = -40 



(24t 3 + 32t 2 + 24t + 5) 3 
448i 7 + 1304t 6 + 1956t 5 + 1735t 4 + 915i 3 + 287t 2 + 51* + 4 



(24t 3 + 32t 2 + 24t + 5) 2 
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Figure 5. Z: T(s) = (s 4 ,s 2 ,s). 



Figure 6. X: y 4 — 2xy 2 + x 2 — xy = 0. 



The curve a(s) = (z 2 {s), z 3 (s)) = (s 2 , s) is a parabola and is not .4(2)- equivalent to X. A curve 
from the family 

/3(6, s) = ( Zl (s) + bz 2 (s), z 3 (s)) = (s 4 + 6s 2 , s) 

is not .4(2)-exceptional and its invariants are given by 



K A \ P = 100 



(36 



Us 2 ) 2 s 2 



56 6 s 2 + 140 s 4 



(6 -14 s 2 ) 3 ' ' AIP (6- 14 s 2 ) 2 

Independently of the value of 6 all curves in the family have the same signature equation 

245 r 3 + 40000 - 448 x 2 + 3780 r x - 1575 r 2 + 14525 x - 6000 r = 0, 

which is different from the implicit equation for the signature for X, and therefore the curves 
from this family are not .A(2)-equivalent to X . We finally consider a two-parametric family 

8(a, s) = (zi(s) + a\ z 3 (s),z 2 + a 2 ^(s)) = (s 4 + a x s, s 2 + o 2 s) 

and find out that for a\ = 20 and a 2 = 2 the implicit equations of the signatures of X and 
the curve parameterized by 5(a,s) are the same. Thus we conclude that Z projects to X by 
a parallel projection over the complex numbers. 

To check that this remains true over M, we look at the invariants of (5(20, 2, s) 



K A \s(20,2,s) 
T A \ s (20,2,s) 



25 (6s 5 + 21s 4 + 84s 3 + 90s 2 + 30s + 25) 2 



(3s 3 + 7s 2 + 25s + 5) 3 



14s 7 + 65s D + 294s° + 535s 4 + 450s d + 415s z + 250s + 25 
(3s 3 + 7s 2 + 25s + 5) 2 ' 



By solving equations K A \s(20, 2, s) = K A \^(t) and ^1^(20, 2, s) = T A \y(t) for s in terms of t, we 
find a real reparameterization s = 1 + At which matches the signature maps of two curves, i.e. 
S\ 7 (t) = S\s(20, 2, 1 + At). Thus the signatures of X and the curve parameterized by 5(20, 2, s) 
are identical and, therefore, Z projects to X by a parallel projection over the real numbers. 

We proceed to find a projection. Since S\j(t) = 515(20, 2, 1 + At), we not only know that the 
exists „4(2)-transformation A that maps j(t) to 5(20,2, (1 + At)), but that for any value of t 
the transformation A maps the point j(t) to the point 5(20, 2, (1 + At). Three pairs of points in 
general position are sufficient to recover an affine transformation. Using three pairs of points 
corresponding to t = 0, 1, 2 we find the affine transformation 



x = 256x + 96y + 21, 



y = lGy + 3 
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that transforms X to the curve parameterized by 5(20, 2, s). From decomposition (9) we can 
now recover a parallel projection 

13 13 113 

x = Z\ Z 2 H Z 3 , V = — Zo H — Z 3 

156 128 32 256 16 8 16 

that maps Z to X. As a side remark, we note that there is no central projection (either complex 
or real) that maps Z to X. 



6 Extensions 

6.1 Projection problem for non-rational algebraic curves 

The projection criteria of Theorems 1 and 2 and the group equivalence criteria of Theorem 3 are 
valid for non-rational algebraic curve. An implementation of the projection algorithm is more 
challenging, however, when the curves are described as zero sets of polynomials, rather than by 
a parameterization. To illustrate these challenges, we can, in fact, consider a rational curve, the 
twisted cubic, but this time we will define this curve by implicit equations. 



Example 9. The twisted cubic is defined as the zero set Z of the system of two polynomials 



gi = z\ — z 2 z 3 and g 2 = z 2 — z\. Following the projection criterion of Theorem 1, we define 



a family of algebraic curves 



Z\ +C\ Z 2 + C 2 

z 3 + c 3 ' z 3 + c 3 



(zi,z 2 ,z 3 ) G Z 



To restrict differential functions, and in particular invariants, to curves from the family i? cliC2iC3 
we need to compute their implicit equations. A naive approach would be to define an ideal 

Y := ((#3 + c 3 )x - (zi + ci), (z3 + c 3 )y - (z 2 + c 2 ) z 3 + c 3 , (23 + c 3 )5 - 1, 
gi^, z 2 , z 3 ), g 2 (z 1 , z 2 , z 3 )) 

and compute the elimination ideal Y n M[ci, c 2 ,c 3 ,x, y]. This leads to a polynomial (cubic in x 
and y) 

A c (x,y) = (c| - ci)y 3 + (c| + c 2 )y 2 x - (c| + c 2 )x 2 + (ci + c 3 c 2 )xy - 3(cic 3 + c 2 3 c 2 )y 2 
+ 2(cic 3 - cl)x + 3(c2C 3 + cic 2 )y -c\- c|, 

which, indeed, defines the curve Z CljC2:C3 provided c| — c\ 7^ and c| + C2 7^ 0, but not otherwise. 
The underlying issue is non-commutativity of specialization of parameters c with the elimination 
(or, in more geometric language, non-commutativity of intersection and Zariski closure). 

To find implicit equations for ,Z CljC2)C3 for the remaining values of the parameters c, we need 
a detailed analysis of the constructible set obtained by the projection of the variety of the ideal Y 
onto (x, y, ci, c 2 , c 3 )-subspace. This could be done using, for instance, RegularChains package 
in Maple. We find out that when c| — c\ 7^ 0, but c 2 + c 2 = 0, the curve Z Cl:C2fi3 is a zero set 
of a cubic 

B c (x, y) = y 3 + 3c 3 y 2 - xy + 3c|y - 2c 3 x + ci + c 3 . 

When both — ci = and c 2 + c 2 = (that is the case when the center of the projection lies 
on the twisted cubic), then Z cltC2;C3 is a zero set of a quadratic polynomial 



C c (x, y) = y 2 + c 3 y -x + c§. 
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On Step 2 of Algorithm 1, where we are describing the family of curves Z CljC2jC3 , we must produce 
all three possible implicit equations A c (x, y) = 0, when c 3 — c\ ^ and c\ + C2 7^ 0, -B c (x, y) = 0, 
when c| — ci 7^ 0, but c 2 + C2 = and C c (x, y) = 0, when c 3 — c\ = and €3 + 02 = 0. Then the 
rest of the algorithm should run for each of these cases with appropriate conditions on c. 

We found that, for majority of the examples, producing the set of all possible implicit equa- 
tions for the curves Z CliC2jC3 (3) from the given implicit equations of an algebraic curve Zcl 3 
to be a very challenging computational task. 

6.2 Projection problem for finite lists of points 

The projection criterion of Theorem 1 adapts to finite lists of points as follows: 

Theorem 5 (central projection criterion for finite lists). A list Z = (z , . . . ,z m ) of m points 
in M 3 with coordinates z l = (z[, z 2 , z\), I = 1, . . . , m, projects onto a list X = (x 1 , . . . , x m ) of to 
points in M. 2 with coordinates x' = (x , y ), I = 1, ...,m, by a central projection if and only if 
there exist ci, 02,03 € M and [A] £ VGC(3), such that 

[x l ,y l ,l] T = [A)[z[ +d,4 + C2.4 + c 3 ] T for l = l,...,m. 

The proof of Theorem 5 is a straightforward adaptation of the proof of Theorem 1. The 
parallel projection criteria for curves, given in Theorem 2 and Corollary 1, are adapted to the 
finite lists in an analogous way. 

The central and the parallel projection problems for lists of m points is therefore reduced to 
a modification of the problems of equivalence of two lists of m points in PR 2 under the action 
of VQC{2>) and A{2) groups, respectively. A separating set of invariants for lists of to points 
in PR 2 under the *4(2)-action consists of ratios of certain areas and is listed, for instance, in 
Theorem 3.5 of [27]. Similarly, a separating set of invariants for lists of m ordered points in 
PR 2 under the 'P<5£(3)-action consists of cross-ratios of certain areas and is listed, for instance, 
in Theorem 3.10 in [27]. In the case of central projections we, therefore, obtain a system of 
polynomial equations on c\, C2 and C3 that have solutions if and only if the given set Z projects 
to the given set X and an analog of Algorithm 1 follows. The parallel projections are treated in 
a similar way. Details of this adaptation appear in the dissertation [6]. 

We note, however, that there are other computationally efficient solution of the projection 
problem for lists of points. In their book [20], Hartley and Zisserman describe algorithms that 
are based on straightforward approach: one writes a system of equations that relates pairs of the 
corresponding points in the lists ZcM 3 and X C M 3 and determines if this system has a solution. 
The book also describes algorithms for finding parameters of the camera that produces an 
optimal (under various criteria) but not exact match between the object and the image. 

In [2, 1], the authors present a solution to the problem of deciding whether or not there exists 
a parallel projection of a list Z = (z , . . . ,z m ) of m points in M 3 to a list X = (x , . . . ,x m ) 
of to points in M 2 , without finding a projection explicitly. They identify the lists Z and X 
with the elements of certain Grassmanian spaces and use Pliiker embedding of Grassmanians 
into projective spaces to explicitly define the algebraic variety that characterizes pairs of sets 
related by a parallel projection. They also define of an object/image distance between lists of 
points ZcK 3 and XcM 2 , such that the distance is zero if and only if there exists a parallel 
projection that maps Z to X. As illustrated by Fig. 7, a solution of the projection problem 
for lists of points does not provide an immediate solution to the discretization of the projection 
problem for curves. Indeed, let Z = (z , . . . , z m ) be a discrete sampling of a spatial curve Z and 
X = (x 1 , . . . , x m ) be a discrete sampling of a planar curve X. It might be impossible to project 
the list Z onto X, even when the curve Z can be projected to the curve X . 
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Figure 7. Projection problem for curves vs. projection problems for lists of points. 



6.3 Applications: challenges and ideas 

A discretization of projection algorithms for curves will pave a road to real-life applications 
and is a topic of our future research. Such algorithms may utilize invariant numerical approxi- 
mations of differential invariants presented in [5, 9]. Differential invariants and their approxi- 
mations are highly sensitive to image perturbations and, therefore, pre-smoothing of the data 
is required to use them. Since affine and projective invariants involve high order derivatives, 
this approach may not be practical. Other types of invariants, such as semi-differential (or 
joint) invariants [27, 32], integral invariants [16, 19, 29] and moment invariants [31, 33] are less 
sensitive to image perturbations and may be employed to solve the group-equivalence problem. 

One of the essential contributions of [1, 2] is the definition of an object/image distance 
between ordered sets of m points in M 3 and M 2 , such that the distance is zero if and only if these 
sets are related by a projection. Since, in practice, we are given only an approximate position 
of points, a "good" object/image distance provides a tool for deciding whether a given set of 
points in M 2 is a good approximation of a projection of a given set of points in M 3 . Defining 
such object/image distance in the case of curves is an important direction of further research. 

A Appendix 

We provide explicit formulae for invariants in terms of jet coordinates. For convenience we recall 
our notation 



Ax :=3 y ( 4 V 2 )-5[ y ( 3 )] 2 , 

A 2 := 9y^[yW] 2 -45y^yWyW+40[yW] 3 . 



(45) 



A classifying set of _4(2)-invariants (23) is given by 




(46) 



(Ai) 3 ' 

92/ (6)[ y (2)] 3 



63y (5) y (3) [ y (2)]2 



45 [y&] 2 [yW] 2 + 255y( 4 ) [y®] V 2 ) - 160 [y^] 4 



(Ai) 2 



while a classifying set of rational 'PC?£(3)-invariants (24) is given by 




(47) 



+ 210y( 5 y 3 )[y( 2 )] 2 



( 



63 [yW] 2 [»«] 2 - 60 9 « [„<»>] V 2 > + 32 [ v 0>] 4 ) 
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- 525 y ( 4 V 2 ) (9 [y^] 3 [y^] 3 + 15 [y^] 2 [y^] 2 [y^] 2 

- 60y( 4 ) [y®] A yW + 64 [y®] 6 ) + 11200 [y^] 8 ]', 

2?/ 8 V 2) (A 2 ) 2 - 8y^A 2 (9y^[y^] 3 - 36y^[y^] 2 



243 [y W 



2(a 2 ; 



4 

45 [y W] 2 [ y ( 2 )] 2 + 12 0y ( 4 ) [y<«>] 2 - 40 [y^] 4 ) + 504 [y( 6 )] 3 [y( 2 )] 5 



504 [y< 6 )] 2 [y( 2 )] 3 ( 9j/< W 2 > + 15 [y^] V 2 ) - 25y ^ [y® 



+ 28y^ (432 [y&] 2 [y( 3 )] 2 [</ 2 )] 3 + 243 [y^] V 4 ) [y^] 4 - 1800y( 5 y 4 ) [y( 3 >] 3 [y^] 2 

- 240y ^ [y^] V 2 ) + 54 0y( 5 ) [y( 4 )] 2 [y( 3 )] [y( 2 )] 3 + 6600 [y^] 2 [y^] V 2 ) 

- 2000 y ( 4 ) [y^] 6 - 5175 [y W] 3 [</ 3 )] 2 [</ 2 )] 2 + 1350 [</ 4 )] 4 [</ 2 )] 3 ) 

- 2835 [y( 5 )] 4 [y( 2 )] 4 + 252 [y^] V 3 ) [y^] 2 ( 9y ( 4 V 2 ) - 136 [y ^] 2 ) 

- 35840 [y( 5 )] 2 [ y ( 3 )] 6 - 630 [y <*>] 2 [ y ( 4 )] [y( 2 )] ( 6 9 [y«] 2 [y^] 2 - 160 [</ 3 )] 4 

- 153yW [y( 3 )] 2 [</ 2 )] ) + 2100y( 5 ) [yW] V 3 ) (?2 [y W] 4 + 63 [„ W] 2 [</ 2 )] 2 (48) 

- 193y W [y^] V 2 >) - 7875 [</ 4 )] 4 ( 8 [y W] 2 [y( 2 >] 2 - 22 y ( 4 ) [y( 3 )] 2 [y ^] + 9 [ y ( 3 ); 
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